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We consider a dark soliton in a Bose-Einstein condensate in the presence of a weak disorder 
potential. Deformation of the soliton shape is analyzed within the Bogoliubov approach and by 
employing an expansion in eigenstates of the Poschl- Teller Hamiltonian. Comparison of the results 
with the numerical simulations indicates that the linear response analysis reveals a good agreement 
even if the strength of disorder is of the order of the chemical potential of the system. In the second 
part of the paper we concentrate on the quantum nature of the dark soliton and demonstrate that 
the soliton may reveal Anderson localization in the presence of disorder. The Anderson localized 
soliton may decay due to quasi-particle excitations induced by the disorder. However, we show that 
the corresponding lifetime is much longer than the condensate lifetime in a typical experiment. 

PACS numbers: 03.75.Lm,72.15.Rn,05.30.Jp 



C/3 



I. INTRODUCTION 

Ultra-cold atomic gases become a playground where 
complex systems of solid state physics, nonlinear quan- 
tum optics or even cosmology can be efficiently simulated 
and investigated [11, Q • The level of experimental control 
and detection is unprecedented and allows one to build 
quantum simulators, i.e. experimentally controlled sys- 
tems that are able to mimic other systems difBcult to 
investigate directly Q . Since the first laboratory achieve- 
ment of Bose-Einstein condensation the list of problems 
investigated experimentally and theoretically in ultra- 
cold atomic gases becomes very long and includes: super- 
fluid phases of both bosonic and fcrmionic atomic species 
11, JH, collective excitations like solitons and vortices @- 
[l2j , insulating phases of solid state physics [l3l - [l^ . tr ans- 
port properties and Anderson localization effects jl6l - [l9l | , 
atomic systems in the presence of artificial gauge poten- 
tials [Utm. 

Particularly interesting is the interplay between par- 
ticle interactions and disorder phenomena in quantum 
many body systems. Anderson localization [23|, which is 
essentially a single particle phenomenon, is vulnerable to 
particle interactions [23 - !29j . In experiments with atomic 
gases the localization was observed when the interactions 
were practically turned off by employing Feshbach reso- 
nances or by reducing density of the atomic cloud [l6l - [l9| . 

Particle interactions play a vital role in soliton forma- 
tion. Mean field description of atomic Bose-Einstein con- 
densates reduces to the Gross-Pitaevskii equation [sol - 
l32j which possesses a localized wave-packet solution (i.e. 
bright soliton) for attractive particle interactions and a 
solution with a phase flip (i.e. dark soliton) for repul- 
sive interactions. Both kinds of the solitons have been 
observed experimentally 

While the Gross-Pitaevskii equation is a single particle 
description with the interactions included in the mean 
held approximation one can anticipate quantum many 
body effects that go beyond such a description. Cen- 
ter of mass of a bright soliton is a degree of freedom 



which, when properly described quantum mechanically, 
allows for the analysis of interesting phenomena like soli- 
ton scattering on a potential barrier which leads to a 
superposition of macroscopically distinct objects [33| or 
quantum entanglement of a pair of solitons |34| . It is also 
shown that in spite of the fact that the particle interac- 
tions are present and are responsible for the bright soliton 
formation the center of mass of the soliton can reveal An- 
derson localization in a weak disorder potential (ssl . [36j . 
One may raise a natural question if similar localization 
can be also observed in the case of a dark soliton. In the 
bright soliton case a huge energy gap for quasi-particle 
excitations guarantees that the shape of the soliton is not 
perturbed by a weak disorder potential. Then, the only 
degree of freedom affected by the disorder is the center 
of mass and (as has been shown recently) it becomes An- 
derson localized. In the dark soliton case the situation 
is more delicate because there is practically no gap for 
quasi-particle excitations. Consequently one may expect 
that coupling between the degree of freedom that de- 
scribes soliton position and the quasi-particle subsystem 
induced by the disorder can destroy Anderson localiza- 
tion of the dark soliton. In the present paper we address 
this problem and show that the dark soliton localization 
can be observed experimentally because the lifetime of 
the localized states is sufficiently long. 

The paper is organized as follows. In Sec. HI] we analyze 
the deformation of a dark soliton solution in the presence 
of a weak external potential within the Bogoliubov ap- 
proach and by reducing the description to the problem of 
the Poschl- Teller potential. Such a classical description 
of the soliton allows us to conclude that the deforma- 
tion effect is very weak. It implies that in the quantum 
description of the soliton we may expect that coupling 
between the degree of freedom that describes soliton po- 
sition and the quasi-particle subsystem can be neglected. 
It is confirmed in Sec. IIIII where we switch to the quan- 
tum description and calculate the lifetimes of Anderson 
localized dark solitons. We conclude in Sec. IIVI 
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II. CLASSICAL DESCRIPTION 

We consider A^o bosonic atoms with repulsive interac- 
tions in a ID box potential of length L at zero tempera- 
ture (for experimental realization of a box potential see 
[37}). The mean field description assumes that all atoms 
occupy the same single particle state 0o which is a solu- 
tion of the Gross-Pitaevskii equation (GPE) [sO] 
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where = 2haLu±, fiQ is chemical potential of the system, 
a stands for the s-wave scattering length of the atoms 
and u}± denotes the transverse harmonic confinement fre- 
quency. By virtue of the nonlinear term in the GPE there 
exists a stationary dark soliton solution which, far from 
the boundary of the box potential, takes the form 



(l)o{x-€j) = e '''^/p^^tanh 
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where 6 is an arbitrary phase, po stands for the atomic 
density away from the soliton position q and ^ = 
h/ ^^mgopo is the so-called healing length. The chemical 
potential of the system is /ig = goPo- The description of a 
dark soliton restricted to the classical wave equation ([T]) 
is called the classical description. It is in contrast to the 
quantum description of Sec. IIIII where, for example, the 
position of a dark soliton q becomes a quantum degree 
of freedom. We consider a finite ID system. However, in 
order to describe the system in a region away from the 
boundary of the box potential we may use analytical so- 
lutions corresponding to the infinite configuration space 
which, if necessary, can be modified to satisfy boundary 
conditions [38|, e.g. 4>o{x — 0) ~ and (j)o{x = L) = 0. 

In the remaining part of the paper we adapt the fol- 
lowing units for energy, length and time, respectively 
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A. Deformation of a dark soliton: Expansion in 
Bogoliubov modes 

We would like to describe the deformation of a dark 
soliton solution in the presence of a weak disorder po- 
tential. The disorder potentials we are interested in cor- 
respond to optical speckle potentials which arc created 
experimentally by shining laser radiation on a so-called 
diffusive plate [H, [s^. In the far field, light forms a 
random intensity pattern which is experienced by atoms 
as an external disorder potential. Diffraction effects are 
responsible for a finite correlation length of the speckle 
potentials. Interestingly properties of such a disorder can 
be easily engineered that allows, e.g., for preparation of 
a matter- wave analog of an optical random laser [40| . 



Bright soliton deformation in an optical speckle po- 
tential has been already considered in Ref. |4l[, see also 
[4^ . Dark soliton propagation and its radiation in the 
presence of randomly distributed Dirac-delta. potentials 
has been considered in Ref. [4^ , see also [13, |4^ . In the 
case of repulsive particle interactions the deformation of 
a condensate ground state in an optical speckle potential 
has been also a subject of scientific publications p6i - |49| . 
In the present paper we consider a problem of the re- 
pulsive interactions but we deal with a localized soliton 
structure, i.e. there is a phase flip of the condensate 
wavefunction at the position of a soliton. 

In our considerations the dark soliton is placed in a 
bounded and weak external potential V{x). To calcu- 
late a small perturbation of the solitonic wavefunction 
we start with the time-independent GPE which in our 
units (|3]) takes the form 
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di(b{x) + -\cj){x)\'Hx) + V{x)cj,{x) = Mx), (4) 
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and substitute 4> = ipQ + 6(/) and p, = p.o + Sfj, = 1 + Sp 
where dcj) is a small perturbation of the soliton wavefunc- 
tion ^ and dp is a small contribution to the chemical po- 
tential which allows us to correct a possible change of the 
total particle number due to the presence of V{x). Keep- 
ing linear terms only we obtain the time-independent, 
non-homogeneous Bogoliubov-dc Genncs equations 
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In order to solve Eqs. ([5]) we would like to expand the 
two component vector ((50, 50*)"^ in a complete basis that 
consists of eigenvectors of the non-Hermitian operator 
C. This basis has been published in Ref. [38|, see also 
[50I . [sH , and here we only present the results and com- 
ment on essential elements of the derivation. Collecting 
the complete basis one has to be careful because the C 
operator is not diagonalizable [H, [H, [13, [EH- There 
are right eigenstates of £, i.e. C\ij)k) = efe|V-'fc)- where 
iV'fe) = \vk)Y and [11] 
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that correspond to the familiar Bogoliubov spectrum 
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These eigenstates are called phonons and their adjoint 
vectors |?/'fc'^) = {\'Uk),—\vk))'^ are also the left eigen- 
modes of C. Due to symmetries of the C operator ~ 
I'Sfe))^ = are the right eigenmodes 

corresponding to eigenvalues — — ej. and — 
{—\vl),\ul.))'^ are their adjoint vectors. In the infinite 
configuration space right eigenvectors and adjoint vec- 
tors fulfin = {uk\uk') - {vk\vk') = S{k - k'). In 

a box potential the wavevector becomes quantized, i.e. 
kn = mr / L where n = 1, 2, . . ., and there is a tiny energy 
gap for the phonon excitations. 

Phonons apart, there exist two zero-eigenvalue modes 
of £ correspondin g to broken gauge U{1) and transla- 
tional symmetries jSq , 
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respectively. They appear as zero-eigenvalue vectors be- 
cause a shift of the global phase 6 or a change of the 
soliton position q in solution ^ does not lead to a 
change of the system energy. The zero modes fulfill 
{uq\us) — {vq\vg) = and they are orthogonal to the 
modes adjoint to the phonons. Vectors adjoint to the 
zero modes are not left eigenvectors of the C and they 
can be found by solving [5l| 
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where Af^.g are determined by normalization conditions 
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where Mg = po{dNo/dp()), 
{2q - L)pQ/MqMg. We have 
and (ug'^'luq) — {vg'^\vq) = 0. A small contribution of 
the zero mode {uq, Vq)^ 
{uf\uf) - {vf\vf) = 



Mq = 



-4po and R = 




in (1131) allows one to fulfill also 



0. 



Now we have all vectors to build a complete basis and 
the deformation of the soliton can be expanded in that 
basis 
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Substituting ([TS]) into ([5]) results in 
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Projecting this equation onto the adjoint vectors we can 
obtain the expansion coefficients and a small correction 
to the chemical potential. Note that, there is no restric- 
tion for a choice of small deviation A6 and Aq which is 
due to the fact that these coefficients are related to the 
zero modes. However, while 9 can be arbitrary, we will 
see that q can not because the external potential breaks 
the translational symmetry but does not affect the U{1) 
symmetry. Coefficient A6 is related to a shift of the 
global phase of the soliton solution ([5]) and without loss 
of generality we may choose AO = 0. For the coefficient 
Pg we get 
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The last term on the r.h.s. can be written in the following 

way 

{Uq\V(^0) + {Vq\V(^l) ^ ( dx\M^~q)\''d,V{x),{lS) 

Jo 

which represents a force acting on the soliton. A similar 
effective force appears in the approach of Ref. [13 . For 
the stationary solution it is obvious that the soliton po- 
sition q is chosen so that such a force is zero. Then, also 
an arbitrary shift of the soliton position should be zero, 
i.e. Aq = in ([T5|) . With such a choice of the soliton 
position the coefficient Pg is directly related to a change 
of the total particle number which we assume to be zero. 
Hence, Pg = dN = and Eq. (fT7|) allows us to obtain 
the correction to the chemical potential 
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Projecting ([T6l) onto the adjoint mode (fM)) results in 
Pq — what can be expected because Pq has an inter- 
pretation of the soliton momentum and for the stationary 
state it should be zero. 

Thus, with a proper choice of the soliton position and 
a suitable correction of the chemical potential, all coeffi- 
cients in (|15p are zero except 
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which contain full information about soliton deformation 
in a weak external potential. Finally the stationary soli- 
tonic solution in the presence of a weak external potential 
reads 



(j>{x) = 00 (x) + ^ [bkUk{x) + blvl{x)] 
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B. Deformation of a dark soliton: Expansion in 
modes of the Poschl— Teller potential 

The Bogohubov approach is suitable in the description 
of the eigenmodes of collective or elementary excitations 
of a Bose- Einstein condensate [3l|, [s^. However, if we 
arc interested in the description of a stationary state of 
the GPE and if such a state can be represented as a real 
function one may introduce a substantial simplification. 
We will see that description of dark soliton deformation 
reduces to an expansion of a wavefunction perturbation 
in modes of the Poschl-Teller potential [5J]- Such an 
approach has been applied in the analysis of a bright 
soliton deformation in a weak disorder potential |4l| and 
it turns out it is also suitable in the dark soliton case. 

Let us start again with the stationary GPE but assume 
the solution we are looking for is a real function 



49l + -0'-M + ^W)0 = O. (22) 
2 Po 



Similarly to the Bogoliubov approach we introduce fi = 
fiQ + Sii ~ l + Sfi and (j) = (j)o + 64> but assume that in ^ 
the global phase 9 = 0. Keeping linear terms only, rewrit- 
ing (l)l{x -q) = Po tanh^(x - q) ^ po [l - cosh"^(x - q)] 
and changing variable x ^ x + q wc obtain 



where 
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is the Hamiltonian for a particle in the Poschl-Teller po- 
tential [5^. To compute 6(j) we need to invert the opera- 
tor Hq + 2. All the eigenstates of the Hamiltonian are 
known in the literature ISJI ■ There are two bound states 



i'oix) = ^sech(x)^. 
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with eigenenergies Eq = —2 and Ei = respectively 
and scattering states 



ipkix) 



k"^ - 2 + 3sech(a;)2 -I- 3ifctanh(a;) 



(27r)i/2 



[(l + fc2)(4 + A:2)]i/2 



(27) 

,2 

with Ek = ^, k G R. We can therefore expand the 
deformation 6(j) over orthonormal basis of eigenfunctions 

Scj) = ao V'o + "1 V-'i + J dk ak ipk{x), (28) 

and compute coefficients aj by projecting Eq. (|23p on the 
proper eigenmodes 



(Ej + 2)aj = / dx V*(x)[<5/i0o - V{x + q)(j)o]. (29) 



The wavefunction ([25]) is a zero mode, i.e. {Ho+2)tfjo = 0. 
Thus, Eq. can be solved provided the projection of 
its r.h.s on the zero mode vanishes. Therefore we require 
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In (PD|) we have taken advantage of (^ol^o) = and the 
fact that the zero mode is also the translation mode of 
the dark soliton, i.e. tpQ ~ dx(l)o- Condition ([5(1)) im- 
plies that the soliton position q should be chosen so that 
the force acting on it is zero, compare with (I18|) . From 
Eq. we do not get any restriction for the value of ao- 
However, because aoV'o has an interpretation of a shift 
of the soliton position we should choose = if we are 
interested in a stationary solution. 

In order to solve Eq. (|23p we have to invert the op- 
erator Hq 2 in the Hilbert space with the zero mode 
excluded which is simple because all eigenfunctions of Hq 
are known. That is 



x) = dy K{x, y) [5p(j)o - V{y + g)0o] 
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where the symmetric kernel K{x, y) reads 
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— 4sh|a; — ylshaish?; — 6|x — 2/1} . (32) 

In the Bogoliubov approach, even if we are restricted 
to the phonon subspace, the C operator possesses a 
zero eigenvalue if the configuration space is infinite, see 
0. Therefore it is not straightforward to obtain a sim- 
ple form of the relevant kernel because integration over 
phonon subspace should actually be substituted by sum- 
mation over discrete values of the wavevector. 

Finally we have to determine the correction to the 
chemical potential. In the Bogoliubov approach there is 
a specific mode responsible for a change of the total par- 
ticle number which is orthogonal to all other Bogoliubov 
modes and to keep the particle number constant it is suf- 
ficient to ensure that the wavefunction perturbation does 
not have any contribution from this mode. In the present 
approach the chemical potential has to be determined by 
the normalization condition (^|^) = A'o + 0{5(t>'^) that 
implies ((/)o|<50) = and thus 
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J dxdy (j)o{x)K{x, y)V{y + q)(f>o{y) 
Jdxdy (j>o{x)K{x,y)(po{y) 
J '^y (tanhy -I- 2/ sech^y) tanhyV^(g -I- y). 
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(33) 

which is the same expression as in the BogoHubov ap- 
proach, compare with (|19p . 

The final expression for a stationary solitonic solution 
in the presence of a weak external potential reads 



4'q{x)- j dy K{x,y)V{y + q)(j)oiy) 
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where we have made use of / dyK{x, y)(f>o{y) — — — 
with (j)Q\fj,o = (j)a being a wave function with fixed chemi- 
cal potential fiQ = I. 

In the bright soliton case [ilj the eigenfunction ipi (not 
Tpo as in the dark soliton problem) turns out to be a trans- 
lation mode of the system and also the Hq + 2 operator 
is substituted by Hq + 1/2. Consequently, one obtains 
a different expression for the symmetric kernel which is 
expected because the bright soliton represents a localized 
wavepacket while a dark soliton is a phase flip at soliton 
position with otherwise uniform density. 



C. Comparison with numerical calculations 

We are going to compare the perturbative approaches 
introduced in the previous subsections with results of nu- 
merical calculations in the case of a dark soliton in a weak 
optical speckle potential. However, first we consider a 
simple harmonic potential 



V{x) = — Vo cos(fcoa;), 



(35) 



which on one hand is a generic example because any po- 
tential can be expanded in a Fourier basis and on the 
other hand allows us to obtain simple expressions for ex- 
pansion coefficients. 

In the Bogoliubov approach Eq. ([20)) . for the potential 
([55)1 . one gets 



2fc(4 + P)Vi^ 
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This coefficient contains Dirac delta functions which tell 
us that the condensate density will be modulated with a 
period of 27r/fco. At small k the coefficient diverges 



l= + oi\k\y'). 



(37) 



but the contribution to the soliton wavefunction from 
long wavelength phonons is finite 




Figure 1: (Color online) In panel (a) we show an example 
of the optical speckle potential with the correlation length 
an = 0.05 and for Vo = 1 while in panel (c) we present the 
corresponding exact solution of the Gross-Pitaevskii equation 
obtained numerically (solid black line) and within the per- 
turbation approach (red dashed line), see Eq. (|34|) (Ea. (|2H) 
leads to the same results). In panels (b) and (d) we show the 
same as in (a) and (c) but for an = 1 and Vo ~ 0.5. Green 
dotted lines in (c) and (d) correspond to unperturbed soliton 
wavefunctions 



Thus, in order to calculate soliton deformation we may 
assume infinite configuration space because the singular- 
ity at fc — > is actually not dangerous. The correction 
to the chemical potential < ^1- (O: disap- 

pears in the limit oi koL ^ oo which is a consequence of 
the self- averaging property of the potential see j47j . 

In the approach that involves modes of the Poschl- 
Teller potential Eq. yields 
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for scattering modes and 
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for the bound state tpi, where we have assumed that the 
system is infinite and thus the correction Sfj, = 0. There 
are Dirac delta functions present in similarly as in 
the Bogoliubov coefficients but contrary to the Bogoli- 
ubov case the coefficients ak possess no singularity at 
k^O. 

Note that the coefficients ([36]) and ([39]) are propor- 
tional to Vo- Thus, for Vo <C 1 (i.e. for the poten- 
tial strength much smaller than the chemical potential 
of the system) the perturbation of the condensate wave- 
function is certainly negligible. Deformation of the soli- 
ton shape depends on the relation of fcg to the inverse 
of the healing length. For fco ^ 1. i.e. if the potential 
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changes on a scale much larger than the heaUng length 
of the system, we may employ the Thomas-Fermi pro- 
file [soj and approximate the condensate wavefunction by 
(l}{x) « a/po[1 - V{x)] tanh(a:/|) where ^ = 1/^1 + Vq is 
a local healing length around the soliton position. Thus, 
for fco ^ 1, the soliton shape is still given by the hy- 
perbolic tangent function but its size is smaller. In the 
other limit, i.e. for fco 3> 1, the condensate wavefunction 
reveals harmonic oscillations superposed on the hyper- 
bolic tangent function. Indeed, for kg ^ 1, Eq. ([5^ im- 
plies that the dominant contribution comes from k w kg 
which behaves like Vb/fcg. Thus, even if we increase Vq 
but choose ko ^ 1 the wavefunction perturbation is still 
negligible due to smoothing effects. Finally we would like 
to stress that plotting <j> obtained in the Bogoliubov ap- 
proach and with the help of the Poschl- Teller modes one 
gets exactly the same results. 

For comparison with numerics we have chosen the case 
of an optical speckle potential. Such a potential is charac- 
terized by: zero mean value V{x) = 0, where the overbar 
denotes an ensemble average over disorder realizations, 

I 1/2 



standard deviation Vo 



Vix) 



and autocorrelation 



function V{x')V{x' + x) = 



2 sin^Jjc/aji) 







(x/aR) 



where is the 



correlation length of the disorder. In Fig. [l]we show ex- 
amples of the solitonic solutions in the presence of an 
optical speckle potential for a correlation length much 
smaller {an = 0.05) and comparable [aji = 1) to the 
healing length and for Vq = 1 and Vq = 0.5 respectively, 
obtained within the perturbation approaches and by a 
numerical solution of the GPE. The agreement is surpris- 
ingly good even though the strength of the disorder is of 
the order of the chemical potential. For an = 0.05 the 
disorder changes rapidly and its effect on the condensate 
is significantly smaller than for aR = 1. 

We have considered a dark soliton in a box potential 
and analyzed its deformation due to the presence of a 
weak disorder potential. Our results can also be applied 
to the system in the presence of, e.g., a shallow harmonic 
trap. Indeed, if we are interested in the deformation of 
the condensate wavefunction in the vicinity of the trap 
center and if the change of the harmonic potential energy 
on a scale of the soliton size is much smaller than the 
disorder strength, i.e. Vq where lj is harmonic trap 

frequency, the effect of the presence of the trap on the 
soliton deformation can be neglected. 



III. QUANTUM DESCRIPTION 

In the previous section, in order to describe ultra-cold 
atoms, we have applied the mean field approximation 
where it is assumed that a many body system is in a 
state where all atoms occupy the same single particle 
wavefunction. Solution of the GPE is an optimal choice 
for such a single particle wavefunction. Then, a station- 
ary dark soliton appears as a solution of the classical 
wave equation and its position is given by a real number 



q. In the present section we will take into account situ- 
ations when particles do not necessary occupy the same 
single particle state. It turns out that the problem can be 
described within the quantum version of the Bogoliubov 
approach where, e.g., q becomes a quantum mechanical 
operator q and the soliton position is described by a prob- 
ability distribution. This approach has rather a semiclas- 
sical nature. The full quantum analysis would involve the 
N-body problem as has been done in Refs. ISSl lSal for the 
bright solitons in optical fibres, see also Ref. j33 |. 



A. Effective Hamiltonian 

An effective Hamiltonian that describes a bright soli- 
ton in the presence of a weak external potential has been 
introduced in Ref. [s^ . It is based on the Dziarmaga idea 
of how to describe non-pcrturbatively de gree s of freedom 
corresponding to Bogoliubov zero modes |38| . Derivation 
of the effective Hamiltonian in the dark soliton case fol- 
lows the same reasoning and therefore we will present key 
elements only. 

In the previous section once the global phase 6 of 
the wavefunction ^ and the soliton position q have 
been chosen no deviations of them have been considered. 
Small deviations can be described with the help of the 
zero modes, compare (Unil-dTT]) and (HI]), while large de- 
viations need modifications of the description. The ex- 
pansion of the wavefunction perturbation around a given 
value of the soliton position q, see is actually not 

necessary because one may treatg as a dynamical vari- 
able and the same is true for 9 [38|. This way we obtain 



(Po 



ad 
,,ad 



,,aa 







+ K 




) 











(41) 



In (|1T|) all modes depend on q and 9 and can follow large 
changes of the soliton position and global phase. Substi- 
tuting (j4ip into the energy functional 



H ^ J dx 
leads to 

H = 



2 2po 



, (42) 



dxV{x)\Mx--q)f' 



p2 

q 

'2\Mg\ ' 
-J2[^kblbk + skibk + bl)], 



with 



Sk 



{uk\VM + {vk\V(l>*o) 



(43) 



(44) 



where only terms of the order 0{P^, b^, PV, bV) are kept. 
Note that Eq. is the Hamiltonian formulation of the 
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classical perturbation theory applied in Sec.IIA. That 
is, fixed points of the Hamilton equations generated by 
(03]) [11] correspond to stationary solutions analyzed in 
Sec.IIA. We know from Sec. II that such stationary solu- 
tions are in a very good agreement with the exact nu- 
merical calculations for the disorder strength even of the 
order of the chemical potential of the system, i.e. for 

In the so-called second quantization formalism the 
quantum many body Hamiltonian corresponds to (|42p 
where the wavefunction (j) is substituted by a bosonic 
field operator (j). Then, also the expansion coefficients in 
(UD) become operators 

Pq - ~idq, (45) 
Pe = N-No = ~idg, (46) 

and fulfill commutation relations 

[q,Pg] = i, 
[0,Pe] = ^, 
[bk,bl,] - Skk'- (47) 

Energy functional (05)) does not depend on 6, thus, in the 
quantum description [Pg , H] = and we may restrict it 
to the Hilbert subspace with exactly Nq particles, i.e. for 
any state in this subspace Pelip) = 0, and the quantum 
effective Hamiltonian reduces to 

H = Hg + HB+Hi, (48) 

where 

Hq = / dxV{x)\M^~q)\' 

^-f^ + ^/d. r(-) V(49) 
Hb = ^Efc^I^fc, (50) 

k 
k 

The Hamiltonian Hq describes soliton motion in an ef- 
fective potential which turns out to be a convolution 
of the original potential with the density |0oP- Own- 
ing to |0oP = /Ootanh^(.T — q) = ^-^^[^ — cosh^(x — q)] 
and V{x) = 0, Hq becomes similar to the corresponding 
Hamiltonian for a bright soliton in a weak external poten- 
tial fssj . The term Hb describes quasi-particle subsys- 
tem (phonons) and Hi is the part of the Hamiltonian that 
couples soliton position degree of freedom with phonons. 
In the classical description (Sec. |ll| such a coupling is 
responsible for the deformation of the stationary conden- 
sate wavefunction. In the following we will not look for 
eigenstates of the total Hamiltonian H but rather con- 
sider eigenstates of Hq and calculate the lifetime of the 



system prepared in these eigenstates due to the coupling 
with the quasi-particle subsystem induced by Hi. 

We would like to emphasize striking differences be- 
tween the classical description and the present quantum 
description. They are most apparent in the absence of 
an external potential. Indeed, for V{x) = 0, the soliton 
position in the classical decription can be chosen arbi- 
trary but it is well defined. In the quantum approach 
the Hamiltonian Hq tells us that eigenstates of the system 
correspond to eigenstates of the momentum operator Pq 
and the corresponding probability distributions for the 
soliton position are totally delocalized. Thus, similarly 
as in the bright soliton case [s^, dramatic condensate 
fragmentations a re p redicted in the quantum approach 
of a dark soliton 15711. 



B. Anderson localization of a dark soliton 

The final form of the Hamiltonian Hq is similar to the 
Hamiltonian for the center of mass of a bright soliton in 
a weak external potential [s^. The effective mass \Mq\ in 
(0^]) is equal to two times the number of particles missing 
in a dark soliton notch while in the bright soliton case 
it is given by a total number of particles in a system. A 
bright soliton is the ground state solution of the GPE 
and excitation of its center of mass increases the energy 
of the system. A dark soliton corresponds to a collec- 
tively excited system and in order to decrease the system 
energy one has to, e.g., accelerate the soliton. Indeed, 
excitations of the soliton position degree of freedom ac- 
tually decrease the system energy due to the minus sign 
in front of the expression (|49p. 

It has been shown that in the presence of a weak disor- 
der potential the center of mass of a bright soliton reveals 
Anderson localization [ssj . The same phenomenon can 
be expected for dark solitons. For V{x) being an optical 
speckle potential with the correlation length smaller 
than the healing length of the system we obtain the ef- 
fective potential in (0^ where the healing length plays 
a role of the effective correlation length. The generic 
properties of the Anderson localization in ID allow us to 
expect that all eigenstates of the Hamiltonian Hq are ex- 
ponentially localized, i.e., have a shape with the overall 
envelope [ssl-lsoj 

l^n(g)l^occxp(^-^^), (52) 

where Hqipn{q) = Enipniq), Qo is the mean position of 
the soliton and koc = hoc{En) is the localization length. 
Indeed, in Fig. [2] we present examples of the Anderson 
localized eigenstates for two values of the standard de- 
viation of the speckle potential Vq. The parameters we 
have chosen correspond to: Nq = 10^ ^^Rb atoms in a 
quasi-lD box potential of L = 3550 (3.37 mm) with the 
harmonic potential of = 2tt x 370 Hz in the transverse 
directions; the correlation length of the speckle potential 
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Figure 2: In top panels we show examples of eigenstates 



panels the corresponding probability densities in log scale. 
The correlation length of the speckle potential aR = 0.28 and 
the strength Vo = 7 x (left panels) and Vo = 1.4 x 10"'' 

(right panels). The eigenstates correspond to the eigenvalue 
En = -3.03 X 10"^ (left panels) and = -8.58 x 10"^ 
(right panels) and reveal the localization length lioc = 10.5 
and lioc = 15.7, respectively. 



CTfl = 0.28 (0.27 //m)[3. The energy units ^ are the fol- 
lowing: Eq/H = 789 Hz, ^0 = 0.96 /im and to ~ 1.27 ms. 

In order to obtain predictions for Anderson localiza- 
tion of solitons we have neglected coupling of the soliton 
position degree of freedom to the quasi-particle subsys- 
tem. In the bright soliton case such an approximation 
is justified because there is a huge energy gap for quasi- 
particle excitations and if the strength of the potential is 
much smaller than the chemical potential of the system, 
corrections to the effective Hamiltonian Hg are negligible 
[36| . In the dark soliton case there is practically no en- 
ergy gap for phonon excitations, i.e. minimal Cfc, see ([9|), 
corresponds to k — n/L which tends to zero for a large 
system. Moreover, a dark soliton is a collectively excited 
state which may decay to lower energy states by emis- 
sion of phonons. If the strength of the disorder potential 
Vo <C 1 we know from the classical analysis that the 
shape of a stationary dark solution of the GPE is negligi- 
bly deformed by the external potential. In the quantum 
description we may thus expect that the lifetime of An- 
derson localized eigenstates is sufficiently long and allows 
for experimental observations of the localization effects. 

Suppose we choose an initial state j^*) of the Aq- 
particle system where the soliton position is described by 
an eigenstate ipniq) of the Hamiltonian Hg corresponding 
to an eigenvalue En and there is no phonon excitation, 
i.e. we deal with the quasi-particle vacuum state. 



I*) 



(53) 



where 6a;|0b) = for each k. In the first order in the 
Hi, see (|5ip . the system may decay to another eigen- 
state ipm (9) corresponding to an eigenvalue Em emitting 
a single phonon of energy e^. According to the Fermi 



golden rule the decay rate reads 

r ^ 27r^7™, 



(54) 



where 



7m = \{i>m,lk\Hi\ipn,0B)\'^9{ek) = \{i>m\sk\ll^n)\g{^k) ■ 

(55) 

The sum in (|54p runs over all eigenstates ipm for which 
we can find such a phonon that the energy conservation 
En = Em+^k is fulfilled. We assume a continuum phonon 
spectrum ^ with the energy gap corresponding to fc = 
tt/L. The density of states is 



[2(e2 + l)(^/^^-l)] 



1/2- 



(56) 



The lifetime of the Anderson localized states presented 
in Figs. Hi and r = l/F = 8 x 10^ (17 minutes) and in 
Figs. [5]d and[2Ji t = 2.5 x 10^ (5 minutes) which means 
that there is by far enough time to perform experiments 
until they can decay due to phonon emissions. In Fig. [3] 
we show contributions jm to the decay rate (j54p and tpm 
states corresponding to the largest values of 7m. Figure[3] 
indicates that the most probable decay leads the system 
to states localized in the vicinity of the initial localization 
region. 

Long lifetime of Anderson localized states is very 
promising from the experimental point of view. Indeed, 
it means that there is sufficient time to excite a dark soli- 
ton in an ultra-cold atomic gas, wait until it localizes in 
the presence of a weak disorder potential and perform 
an atom density measurement. If the soliton is Ander- 
son localized distribution of soliton positions collected in 
many realizations of the experiment [57| will reveal an 
exponential profile. 

Experiments with dark solitons have been performed 
in the presence of a harmonic trap 0, l60l - [6^ . In such a 
case, in order to observe the Anderson localization of the 
solitons, the trap has to be sufficiently shallow. That is, 
the ground state extension of the soliton position in the 
harmonic trap without the disorder must be much larger 
than the localization length predicted in our analysis, i.e. 

, ^ ^ hoc where w is the harmonic trap frequency. 



IV. CONCLUSIONS 

We have considered a dark soliton in dilute ultra-cold 
atomic gases in the presence of a weak disorder potential. 
Our consideration is divided into classical and quantum 
descriptions. The classical approach concerns the anal- 
ysis of stationary solutions of the Gross-Pitaevskii equa- 
tion and the effect of deformation of the soliton shape 
by the disorder. We have employed two methods: the 
Bogoliubov approach and the expansion of a wavefunc- 
tion perturbation in eigenmodes of the Poschl- Teller po- 
tential. These two methods lead to the same results, 
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Figure 3: (Color online) Top panels: contributions 7m, 
Eq. (|55p . to the total decay rate F as a functions of energy 
Em- Bottom panels: initial states i^niq) of the system (solid 
black lines) and eigenstates il^miq) (dashed red lines) corre- 
sponding to the most probable decay channels, i.e. maximal 
values of 7m- Parameters chosen in left (right) panels are the 
same as in the corresponding panels in Fig. [21 



however, the expansion in the Poschl- Teller modes turns 
out to be more convenient and in particular allows us 
to obtain a very simple form of the soliton perturbation 
in terms of an integral kernel. Comparison of the per- 
turbative calculations with the numerical results shows 
surprisingly good agreement even for the strength of an 



external potential as great as the chemical potential of 
the system. If the strength is much smaller than the 
chemical potential the wavcfunction deformation is neg- 
ligibly small. 

The Bogoliubov approach is invaluable in the quan- 
tum description where we are interested in many body 
eigenstates of the system. If the strength of an exter- 
nal potential is much smaller than the chemical potential 
the dark soliton position may be described by an effec- 
tive quantum Hamiltonian which is weakly coupled to 
the quasi-particle subsystem. The effective Hamiltonian 
turns out to be similar to the corresponding Hamilto- 
nian in the problem of a bright soliton in a weak external 
potential. Similarly as in the bright soliton case we pre- 
dict Anderson localization of a dark soliton in the pres- 
ence of a disorder potential. Because there is a coupling 
between the soliton position degree of freedom and the 
quasi-particle subsystem the localized states may decay 
via phonon emission process. We have investigated life- 
times of the Anderson localized states and it turns out 
that for typical experimental conditions they exceed con- 
densate lifetimes that make experimental observations of 
the dark soliton localization realistic. 
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